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Abstract 

In this work we report preliminary results on the relaxational dynamics of one dimensional Bose 
described by the Lieb-Liniger model, upon release from a parabolic trap. We explore 
the effects of integrability and integrability breaking upon these dynamics by placing the gas post- 
release in an integrability breaking one-body cosine potential of variable amplitude. By studying 
the post-quench evolution of the conserved charges that would exist in the purely integrable limit, 
we begin to quantify the effects of the weak breaking of integrability on the long time thermalization 
of the gas. 



I. INTRODUCTION 



In their seminal experiment on one dimensional interacting cold atomic gases, Kinoshita 
et al. [I] argued for the possibility that the relaxational dynamics of such gases possessed 
memory of the gases' initial condition. Specifically, they observed that the momentum 
distribution of the gas did not rapidly evolve to a thermal equilibrium state, despite the 
presence of interactions between the gas' atoms. To explain this behavior, they conjectured 
that the gas possessed non-trivial conserved integrals of motion (beyond the energy of the 
gas), and that these integrals of motion were controlling the long time dynamics of the 
gas. These non-trivial integrals of motion should exist, they argued, because the underlying 
theoretical description of the gas, the Lieb-Liniger model, is known to be exactly solvable and 
has an infinite set of such integrals [21 [3]. In subsequent work Rigol et al. [3] sharpened this 
conjecture by arguing that while the gas did relax to a state governed by a thermodynamic 
ensemble, this ensemble was not the canonical (or microcanonical), but an ensemble aware 
of these additional integrals of motion, an ensemble they dubbed the generalized Gibbs 
ensemble. 

In subsequent work this ensemble has been shown to govern the dynamics of a number of 
systems characterized by sets of non-trivial conserved quantities, both non-interacting |5]-|9] 
and interacting |10l [11] . However less studied has been the question of thermalization when 
the system has a set of weakly broken integrals of motion [T2HI1]. Does the weak breaking 
of the integrals of motion always lead to eventual thermalization of the gas as governed by 
the canonical ensemble? Is there a time scale of integrability breaking, tjb, for which for 
times t < TjB the dynamics appear integrable while for times t > tjb, the dynamics are 
governed by the standard thermodynamic ensembles? Or is there a smooth crossover in 
thermalization as suggested in [T2] - [Ti] . where physical quantities interpolate between their 
values under the generalized Gibbs ensemble and their values in the canonical ensemble? 

It is with these questions in mind that we study the following quantum quench problem. 
We begin by considering a one dimensional atomic Bose gas of N particles in a system of 
length L in the presence of a one-body parabolic trap, Vpa.ra{x) = We describe 

the gas with the Hamiltonian, 
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where we work in units of 2m = h = 1. Hn is the well known Lieb-Liniger model [2] with 
interaction strength c. We prepare the system in the ground state of this Hamiltonian, H. 
At t = we remove the parabolic trap. In the absence of the parabolic trap, the Hamiltonian 
(now just the Lieb-Liniger model itself) is integrable, and we expect the subsequent dynamics 
of the gas to be governed by the non-trivial conserved charges in the system. If instead 
of simply removing the trap, we replace it with a different one-body potential, V^os(a^) = 
Acos(co'x), we break integrability, so changing the nature of the post-quench dynamics. By 
varying the amplitude A of this potential, we control the amount of integrability breaking in 
the system and its concomitant effects on the dynamics. This quantum quench is illustrated 
in Fig. 1. 

To study these dynamics, we employ a combination of Bethe ansatz solvability and a 
numerical renormalization group. The Lieb-Liniger model is a model that can be solved with 
the Bethe ansatz, both determining its spectra [2J and its matrix elements [15]. Because we 
can compute matrix elements, we can compute correlations functions of this model using 
their Lehmann representations [121 HZ] . This however is computationally intensive. In order 
to accomplish this task we employ an optimized set of routines known as ABACUS [16]. 

Because we are interested in perturbing the gas by introducing integrability breaking 
one-body potentials, both pre- and post-quench, the solvability of the Lieb-Liniger model is 
insufficient for the task at hand. Instead to study such deformations, we employ a numerical 
renormalization group (NRG) [TH [THj-EI] able to study perturbations of integrable and 
conformal continuum field theories. This approach, as it is an extension of a methodology 
known as the truncated conformal spectrum approach [221 [23] . has been primarily used to 
study perturbations of relativistic field theories [T8I42T] . but has recently been applied to 
the Lieb-Liniger model perturbed by a one-body potential [11], the problem at hand. The 
NRG uses the eigenstates of the Lieb-Liniger model as a computational basis. Because this 
basis accounts for the interactions of the Bose gas particles with one another, this numerical 
method builds in the strong correlations present in the problem right at the start. 

The paper is organized as follows. In Section 2 we demonstrate that we can compute 
the equilibrium properties of the gas in the one-body potentials. While we have shown in 
Ref. [11] that we can accurately describe the ground state and first few excited states in 
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FIG. 1: Sketch of the quench protocol. For t < 0, we prepare the gas in its ground state in a 
parabolic potential. At t = 0, we remove the parabolic potential and replace it with a small cosine 
potential. 

such a potential, here we show that we can obtain with reasonably accuracy a wide range 
of the spectrum. This will be important for the determination of the post-quench time 
evolution of the Lieb-Liniger conserved charges. In Section 3, we consider the evolution 
of the expectation values of these charges. This will allow us to begin to understand the 
consequences of integrability breaking. In Section 4, we discuss these results briefly and 
examine possible further directions for this work. 

II. EQUILIBRIUM PROPERTIES OF THE GAS IN THE ONE-BODY POTEN- 
TIALS 

In this section we demonstrate that we can compute the equilibrium properties of the gas 
in the one-body potentials necessary to describe the dynamics of the system post-quench. 
We first show that we can prepare the initial state, the ground state of the gas in the presence 
of the parabolic potential, V^ara? accurately. To this end we show in the central panel of 
Fig. 2 the density profile, p(x), of the gas in this trap. We see that we get good agreement 
between the NRG numerics for = L = 14, c = 7200 and an analytical computation of the 
density profile of the gas in its Tonks-Girardeau limit (c = oo). Such analytics are possible 
because the gas, for certain quantities such as p{x), can be treated as equivalent to free 
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FIG. 2: Right: The parabohc one-body potential, Vpa_ra{x) = 0.9x^ (red), used to prepare the 
initial state and the post-quench one-body cosine potential, Vcos{x) = cos{4ttx/ L) with amplitude 
A = 1 (black). Center: The density, p{x), of the gas {N = L = 14, c = 7200) in its ground state 
in the presence of the parabolic potential: red (analytics), black (NRG). Right: The density, p{x), 
of the gas {N = L = 14, c = 7200) in its ground state in the presence of the cosine potential with 
amplitude A = 2: red (analytics), black (NRG). 

fermions. In Fig. 2 (right panel) we also demonstrate that we can accurately compute the 
density profile of the gas in a cosine potential. 

However in order to compute post-quench dynamics, we need to be able to describe not 
only the ground state in the cosine potential, but some large number of excited states. In 
our quench protocol, we take as our initial t = state the ground state of the gas in the 
parabolic potential, \ipGS, para.)- If we can compute a wide range of eigenstates in the cosine 
potential, both ground and excited states, |V'«,cos), we can expand this initial state in terms 
of the post-quench basis: 

(2) 

a 

Of course for this expansion to be exact, we would need to know all of the eigenstates of the 
gas in the parabolic potential. We will instead settle for a determination of the post-quench 
eigenbasis that allows us to include enough states so that |cq,P > 0.99. 

Because we use the eigenstates of the Lieb-Liniger model absent a one-body potential, 
\4'a,LL), as the computational basis of the NRG, the NRG gives any eigenstate in a one-body 
potential as a linear combination of such states: 

I ^one— body 

a 

Each Lieb-Liniger state \ip)a,LL is characterized by A^-rapidities, Aj, i = 1, . . . , N. These ra- 
pidities govern the energy, E^, and momentum, P^, of the state (relative to the Hamiltonian, 
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FIG. 3: A plot of the energies of the gas {N = L = 14, c = 7200) in the cosine potential as 
determined by the NRG (black) and by analytics (red) in the hardcore limit (c = oo). 
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These rapidities are found as solutions of the Bethe equations: 
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In the limit of c = 0, we see that the Bethe equations collapse to the momentum quantization 
condition for a particle in a periodic system of length L. 

In computing the spectrum of states in the cosine potential, we employ the variant of the 
NRG discussed in Ref. pO]. The NRG in its plain vanilla formulation [18] can compute the 
spectrum of the low lying states of the gas in the one-body potential [H]. But to capture 
accurately an appreciable fraction of the spectrum, we need to employ a sweeping routine 
[20] analogous to that used in the finite volume routine of the density matrix renormalization 
group [21]. 

In Fig. 3 we present the results for the spectrum of the gas in the Kos(a^) = cos(47rx/L) 
[A = 1) as computed with the NRG and with analytics in the c = 00 limit. We see we are 
able to describe accurately a wide range of the spectrum. For higher energy states, there 
are some slight differences between analytics and the NRG which we believe can be ascribed 
to 1/c corrections, which while small are still present. 
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FIG. 4: The time evolution of the conserved charges, of the Lieb-Liniger model post-quench. 
Shown are Q/^^ Qq, and Qio for a release into a cosine potential of amplitude A = 0.1,2, and 5. 
The frequency of the cosine potential, u, is set to An /L. 

III. POST-QUENCH DYNAMICS OF CONSERVED CHARGES 

In this section we consider the time evolution of the Lieb-Liniger conserved quantities, 
Qn- These quantities commute with the Lieb-Liniger Hamiltonian, [Qm-f^ii], but in the 
presence of the one-body cosine potential they become time dependent. To compute this 
time evolution, we first note the time evolution of our initial t = state is expressible as 



(6) 



To compute 

Quit) = (V'G5,para(i)|Qn|V'GS,paraW), 

we then need to know (V'/3,cos|Qn|^a,cos)- But because iV'o.cos) is given in terms of Lieb-Liniger 
states (Eqn. |3|, this amounts to knowing the action of the charges, Qn, on such states. This 
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FIG. 5: Amplitudes of the oscillations of the conserved charges, Q4, Qq, and Qio, as a function of 
the amplitude, A, of the post-quench cosine potential. 



however is straightforward j3]: 
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The first two charges in the sequence, Qi and Q2, give the momentum and energy respectively 
of the Lieb-Liniger state. Because of the system's parity symmetry, x — )■ —x, all of the odd 
charges evaluate to zero on \ipLL,a)- We thus will focus on the even charges, Q2n, alone. 

In Fig. 4 we present the time evolution of three charges Q4, Qq, and Qio for quenching 
into cosine potentials with amplitudes A = 0.1, 2, and 5. We see the charges oscillate in time 
with increasing amplitude as the amplitude of the integrability breaking cosine potential is 
increased. These oscillations occur about a well defined mean. This mean smoothly increases 
from its A = value as we increase the strength of the cosine potential. We also see for 
small times, the expectation values of the charges are characterized by transients, but then 
settle into a steady state oscillatory behavior. 

In Fig. 5 we plot the amplitudes of the oscillations of the charg function of the 

amplitude of the cosine potential. We see that these amplitudes of the oscillations are 
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FIG. 6: Average frequency of the oscillations of the conserved charges, Q4, Qq, and Qiq as a 
function of the amplitude, A, of the post-quench cosine potential. 

(roughly) linearly related to the amplitudes, A^s, of the cosine. In the limit that A vanishes, 
the charges become constants of motion, as expected. More interestingly, however, the 
frequencies at which the charges are oscillating. Fig. 6, are independent of the amplitude of 
the cosine potential. Instead these frequencies are roughly but uniformly equal to that of 
the cosine potential itself, in/L ~ 1. 



IV. DISCUSSION AND FUTURE DIRECTIONS 



In this brief report we have considered post-quench dynamics of a cold atomic gas 
quenched from its ground state in a parabolic trap to a cosine shaped trap. To character- 
ize these dynamics we have investigated the time dependency of the Lieb-Liniger conserved 
quantities induced by the presence of the post-quench integrability breaking cosine. In the 
presence of integrability breaking, the charges oscillate about a mean. This mean behaves 
smoothly as a function of the strength of integrability breaking, the amplitude. A, of the 
cosine. Similarly the amplitude of the oscillations go smoothly to zero as A is reduced to 
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zero. However for the frequency of the oscillations of the charges this is not the case. In- 
stead this frequency is related directly to the cosine potential's own frequency, which in our 
quench protocol is kept constant. 

Our findings for the behavior of the charges are then in accordance with Ref. [T2] - [Ti] . 
namely we find that the introduction of integrability breaking leads to a smooth interpolation 
of the expectations values of observables (here the charges) between that in the integrable 
limit and that in the fully chaotic limit where the standard thermodynamic ensembles govern 
dynamics. However we also find that the time scale, r/^, for integrability breaking is not 
necessarily associated with the strength of the integrability breaking. We find rather that 
the charges oscillate with a frequency independent of the amplitude of the one-body cosine 
potential. This suggests the intriguing possibility that even relatively large integrability 
breaking terms, if suitably low-frequency, might lead to long thermalization times before 
completely chaotic behavior is observed. 

In future work we intend to explore these questions in terms of the momentum distri- 
bution function (MDF) of the gas. The long time behavior of the MDF is considerably 
more complicated to compute than the time dependency of the charges. In the presence of 
integrability breaking, it requires one to compute a large number of matrix elements of the 
form, 

(a,LL|V'(a;)V>t(0)|/3,LL), 

where tplx) is the Bose field operator. Each of these matrix elements is evaluated by inserting 
a resolution of the identity between the fields ip, itself a computationally intensive task [T7] . 
Having already carried out preliminary computations for systems sizes of = L = 14, 
we nonetheless expect to be able to perform these calculations for system sizes of up to 
~ 25 - 30. 
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